Spontaneous generation of spin-orbit coupling in magnetic dipolar Fermi gases 
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The stability of an unpolarized two-component dipolar Fermi gas is studied within mean-field 
theory. Besides the known instability towards spontaneous magnetization with Fermi sphere de- 
formation, another instability towards spontaneous formation of a spin-orbit coupled phase with 
a Rashba-like spin texture is found. A phase diagram is presented and consequences are briefly 
discussed. 



Recently, Bose-Einstein condensation of Dysprosium 
164 has been achieved [H. ^^^Dy is an atom with mag- 
netic dipole moment lO/is (/is being the Bohr magne- 
ton). Also a dipolar Fermi gas was recently produced 
in the form of a gas of fermionic molecules with large 
electric dipole moment This gas was, however, not 
degenerate. But since there exists a fermionic Dy isotope, 
i63j3y^ we can expect in the near future the realization 
of a degenerate Fermi gas of atoms with large magnetic 
moment. 

Dipolar atoms interact via the dipole-dipole interac- 
tion (DDI), which is anisotropic and long-ranged. In 
Fermi gases, this interaction has the effect to deform the 
Fermi surface. This phenomenon has been studied in 
Refs. In nuclear physics, the so-called tensor force 

has a similar structure as the DDI Q . The effect of the 
tensor force in nuclear matter has been investigated in 
the framework of Landau Fermi liquid theory, e.g., with 
regard to the spin susceptibility 0, Q . 

Since their recent realization [9], cold atomic systems 
with artificial spin-orbit coupling (SOC), in two as well 
as three dimensions, have received tremendous attention 
(see, e.g., l3"l3|)- In this letter, we show that in three 
dimensional unpolarized dipolar Fermi gases, the DDI 
can give rise to an instability towards spontaneous for- 
mation of a phase with SOC. This could be an alternative 
way, maybe simpler than the artificial SOC, to produce 
SOC in ultracold Fermi gases, and opens wide perspec- 
tives which have been intensely discussed in the very re- 
cent literature. For instance, SOC may have important 
consequences for pairing [ll. 12|. 

The possibility of a spontaneous generation of a SOC 
phase will be studied within the mean-field and random- 
phase approximation (RPA) approach. We will first for- 
mulate the theory for a magnetic dipolar Fermi gas with 
arbitrary spin s and then specialize to s = 1/2. 

The magnetic DDI between two atoms with dipole mo- 
ments di and d2 is given by 
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Let us consider a uniform gas of fermionic atoms hav- 
ing a magnetic dipole moment d — daS, where s is the 
spin operator and do characterizes the magnitude of the 
dipole moment. The Hamiltonian is given by (we use 

units with h = I) 
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The interaction includes the contact interaction as well 
as the dipole-dipole interaction: Ki^,Q'^'(q) = V^p ^,p, + 

The contact interaction is written as Vj^^ ^/^, = 
gSaa'Spp'- For simplicity we assume it to be spin in- 
dependent, although in general its strength depends on 
the total spin of the two interacting atoms [l^ . This ap- 
proximation becomes exact in the case of s = 1/2 atoms, 
which we shall discuss below, since in this case only the 
spin-singlet channel contributes. 

The dipole-dipole interaction V'^'^{q) is the Fourier 
transform of Eq. ([1]), 
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where Cdd = t^o denotes the coupling constant of the 
dipole-dipole interaction and Saa' is the matrix element 
of the spin operator between the basis spin functions 

{Sm)aa' = \/s(s -|- 1)C|°,;^„ with S±i = T{Sx ± iSy) / V^, 
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and the Clebsch-Gordan coefficient C|",j„j 



where r is the distance between the atoms. 



the notation of Ref. [ll| . 

We assume for the moment that the ground state of the 
system is spin symmetric. In this case, the DDI does not 
contribute to the mean field, which implies that all spin 
components have the same spherical Fermi surface (note 
that, as shown in 0, the Fermi surface deforms if the 
ground state is spin asymmetric, i.e., if the gas is fully or 
partially polarized). Then the occupation numbers are 
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given by Pa/3(k) = (4_^Ck,a) = Sapp{k) = SapO^kp - 
k), where kp denotes the Fermi momentum. The single- 
particle energies are — k"^ /(2m) + 2sgn, where 2sgn 
is the Hartree-Fock mean field, n = kp /{6Tr'^) being the 
density per spin state. 

We will use RPA theory in order to investigate for 
which parameters the symmetric ground state gets un- 
stable with respect to the formation of more interesting 
asymmetric phases. To that end, we calculate the spec- 
trum of the collective zero-sound modes: A vanishing or 
imaginary frequency indicates an instability. Let us con- 
sider the retarded response function 

tU^l^^a, ^, {k,k' ,q,t - t') = 0{t - t') 

X ([cL^(i)ck+q.c.W,4,+q_„,(t')ck',/3'(i')]>- (4) 

Within RPA, it is obtained as the solution of the integral 
equation 

n„^,„.;3'(k,k',q,Cj) = nO(k,q,L^)(2^)3,5(k-k')(5aa'5/3/3' 
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X nai^;i,a'/3'(ki,k',q, , (5) 



where /q/3,q'/3' (q, k - k') = Vai3,a'l3'i<i) - Vap,p'a'{k - 

k') denotes the antisymmetrized matrix element of the 
interactionandn*'(k, q, a;) ~ [p(k)— p(k+q)]/(tj— eq+k+ 
£k + iv) is the non-interacting response function. 

In the limit g ^ kp^ the response function is con- 
centrated a,i k = k' ~ kp and it depends only on the 
directions of k and k' and on the dimensionless quantity 
w = muj / {kpq) (without loss of generality we suppose 
that q — qez). As usual in Landau Fermi liquid theory, 
we expand the angular dependence in spherical harmon- 
ics Ylml ■ A complication arises from the fact that the 
DDI is not diagonal in orbital angular momentum L and 
spin S of the excitation, but only in the total angular mo- 
mentum J. We introduce a multi-index A = LSJM and 
define the angular momentum projected response func- 
tion as 
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Then the RPA equation ([5]) reduces to a matrix equation: 

HaaK^^) -II^A'(^)+ E nAA,('^)i^A,A,nA,A'(^), (7) 
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where 
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The Landau parameters -Faa' can be decomposed into 
the direct [D) and exchange [Ex) contributions from the 
contact (c) and dipole-dipole [dd) interactions, i<AA' = 
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explicit expressions read: 
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where the standard 3j, 6j and 9j symbols |14| have been 
used, and Hq = 0, iJ„ = J2p=i ^/P harmonic 
numbers. 

At this point it seems necessary to add a comment 
about the direct term of the DDI in Eq. dTTI) . As one 
can see from Eqs. ([9])— (|12p . the contribution of the di- 
rect term of the DDI is the only one that depends on 
M. Since the Landau parameters are defined in the 
limit q — >■ 0, one would expect that there should not 
be any preferred direction and therefore no dependence 
on M. In fact, for q = 0, one sees immediately that 
ydd^^ _ _ J (fS^y^^^j.^ _ g^j-^^ therefore the direct 

term was omitted in the analysis of Ref. [15]. However, 
the DDI is discontinuous at q = because of its long- 
range nature, and as one can see from Eq. ([3]), V'^'^{q) 
depends only on the direction of q and therefore does 
not vanish in the limit q — >■ 0. Since an instability occurs 
when the energy of a zero-sound mode (which exists only 
at small but non-zero q) vanishes, we will include the 
direct term, as it was done in Ref. |16j . 

Equation ([7]) can now easily be solved by matrix inver- 
sion: 
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where 



A^AA'(t^) = ^AA' -'^'[^AA^i'^)PA^^' 



(14) 
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The onset of an instabihty is characterized by a vanish- 
ing excitation energy. We therefore consider the response 
function in the case co = 0. In this case, Eq. ^ reduces to 
nAA'(O) = -'5aa', so that 7Waa'(0) = Saa' + Faa'- Sta- 
bihty requires that aU eigenvalues of A^(0) are positive 
[fi, 16., J/D| . Note that the Landau parameters Faa' and 
hence the stabihty matrix Ai are diagonal with respect 
to the total angular momentum J. 

Let us now discuss the stability conditions for the spe- 
cial case of s = 1/2 atoms. In this case, the total spin S 
of the excitations can be 5* = or 5 = 1. 

In the channel J — 0, there exist two uncoupled modes: 
A = LSJM = 0000 and 1100. The A = 0000 mode is 
associated with the compressibility and an instability in 
this channel means that the system will collapse. The 
corresponding stability condition is independent of the 
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and reads 
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For A = 1100, the only contribution comes from the 
exchange term of the DDI, since i^noo\ioo — -^iito^iioo — 
-P'noMioo = 0- From Eq. ^ one obtains F^^^^JI^q = 
inkpCdd/i'^T^), and since Cdd is positive, this mode is al- 
ways stable. 

In the case of J = 1, there exist four modes: A — 
lOlM, lllM, OllAf and 211M, the last two being cou- 
pled to each other. For A = lOlM, no instability can 
occur because Fioio.ioio = 0. 

The response in the L = 0, 5* = 1 (i.e., A = OllAf) 
channel is related to the spin susceptibility 0, 13] , and an 
instability in this channel therefore is an instability to- 
wards spontaneous magnetization. Since the A = Oil A/ 
and 211 A/ channels are coupled, the stability condition 
is that the matrix 
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has positive eigenvalues. This condition depends on 
both g and Cdd- The relevant non- vanishing contribu- 
tions from the exchange term of the DDI [Eq. (fT2|) ] are 

Foiim%im/^ = F211M.I11M = rnkpcdd/ {I27r). Note 
that the stability condition depends on M because of 
the direct term of the DDI [Eq. (fTTj) ]. Physically, the 
M = ±1 and M = modes correspond to transverse 
and longitudinal spin waves, respectively. For M = ±1, 



F 



dd(D) 



OHM OHM is negative, while for M = it is positive. 
Therefore the M — ±1 modes become unstable before 



an instability appears in the q 
rect term of the DDI. 



case without the di- 




FIG. 1: Schematic representation of the spin texture induced 
by the perturbation (|18p or created spontaneously if the sys- 
tem is unstable in the A = lllM channel. The circle repre- 
sents a section though the Fermi sphere. The vector a points 
towards the reader. 



For A = lllM, the only non- vanishing Landau pa- 



rameter is F\ 
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111A/,111M 



which is independent of M . The 



stability condition for this channel reads 
rnkpCdd < 47r. 



(17) 



Let us discuss the physical meaning of this instability. 
The mode A = lllM can be excited by the following 
perturbation 



/d'^k 
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where a characterizes the amplitude and orientation of 
the perturbation. This perturbation tends to turn the 
spin s of an atom into the direction perpendicular to the 
momentum k and to the vector a. As a consequence, 
the Fermi surfaces of atoms whose spins point into the 
favored and unfavored directions split up and a spin tex- 
ture in momentum space as shown in Fig. [1] is created. 
This is very similar to the so- called Rashba SOC 
which has recently been discussed also in the context of 
(two dimensional) systems of ultracold atoms [l3|- If an 
instability occurs in this channel, this means that even 
in the absence of the perturbation ^TE\\ . the system will 
spontaneously choose a direction a and create a spin tex- 
ture with spins perpendicular to a and tangential to the 
Fermi surface, as shown in Fig. [T] 

The different stability conditions are summarized in 
Fig. [5] which shows the stable regions as functions of the 
coupling constants g and Cdd- The long-dashed vertical 
line indicates the boundary satisfying Eq. (fT5|) : in the 
region to the left of this line, the mean- field theory pre- 
dicts a collapse. The short-dashed curve indicates the 
onset of the instability towards magnetization obtained 
from Eq. (jl6[) . Finally, the region above the horizon- 
tal solid line, obtained from Eq. p7)) . is unstable with 
respect to the spontaneous formation of SOC discussed 
above. Therefore, the stable region of the spin symmet- 
ric ground state is the gray area. The region where one 
can expect to find the spontaneous formation of SOC is 
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FIG. 2: The phase diagram for the stabihty of the spin sym- 
metric ground state of the system as function of the couphng 
constants g and Cdd- The hues indicate the onset of different 
instabihties: coUapse (long-dashed vertical line), spontaneous 
magnetisation (short dashes), and spontaneous formation of 
a Rashba-like spin texture (solid horizontal line). The spin 
symmetric ground state is stable in the gray area. 



situated above the thick part of the solid line. For com- 
pleteness, it should be mentioned that the modes with 
J > 2 do not change the phase diagram: For J > 2, 
the contact interaction does not contribute and the crit- 
ical values for Cdd are always larger than that given by 
Eq. (dZD. 

Notice that the present stability analysis in the frame- 
work of RPA is only able to detect second-order phase 
transitions. To investigate the stability against first- 
order phase transition, one would have to use other meth- 
ods, e.g., a variational ansatz as in Ref. [4]. 

So far we have only considered the response at uj — 0. 
But using the RPA, we can also calculate the energy spec- 
trum of the corresponding zero-sound modes. This is 
somewhat more difficult because at cj 7^ the response 
function n^A' is no longer diagonal with respect to J 
(although it remains diagonal with respect to M) and, 
therefore, an infinite number of modes are coupled among 
one another in the matrix A4{uj) of Eq. To calculate 
numerically the response function, it is necessary to trun- 
cate the matrix at some value Jmax- However, we found 
that convergence is practically reached at Jmax ~ 10. As 
an example, we display in Fig. |3] the imaginary part of 
the response function in the channel A = 1110 for differ- 
ent values of the coupling constant Cdd- Since the Landau 
parameter -Fmo.iiio is always negative, there is no zero- 
sound mode but only a broad particlc-holc continuum. 
As one approaches the instability at mkpCdd — 47r, the 
response gets more and more enhanced at low energy and 
at the instability it diverges at a) = 0. 

To summarize, we have applied mean- field and RPA 
theory to an unpolarized dipolar Fermi gas in the spin 
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FIG. 3: The imaginary part of the response function (|13|l for 
A = A' = 1110 for various values of Cdd- 

symmetric ground state. We have discussed the stabil- 
ity conditions in the special case of spin-1/2 atoms. In 
addition to the known collapse and spontaneous magne- 
tization instabilities, we found that for certain values of 
g and Cdd the system gets unstable towards a phase with 
SOC where a Rashba-like spin texture around a sponta- 
neously chosen axis in momentum space is formed. 

An important subject for future studies will be to check 
whether the SOC phase survives also beyond the mean- 
field approximation. For instance, in reality the collapse 
predicted by the Hartree-Fock approach in the case of a 
strongly attractive contact interaction does never occur 
and the Fermi gas stays stable even in the unitary limit 
[20I . The stabilizing effect of short-range correlations was 
also discussed in the context of the instability towards 
spontaneous magnetization in the case of a strongly re- 
pulsive contact interaction (without DDI) [2l[ . In nuclear 
physics, short-range correlations due to the tensor force, 
which is similar to the DDI, are known to be important 

In addition, it is interesting to see what happens if the 
gas is trapped and not uniform. It will also be important 
to study not only the onset of the instability, but also the 
SOC phase itself and the competition between the differ- 
ent phases. Another interesting question is whether the 
SOC phase persists in (quasi-) two dimensional systems. 
In this case it would be possible to study, e.g., the spin 
Hall effect with cold atoms. 

We thank D. Basko for interesting discussions and sug- 
gestions. 
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